Decorrelating a compressible turbulent flow: an experiment 
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Floating particles that are initially distributed uniformly on the surface of a turbulent fluid, 
subsequently coagulate, until finally a steady state is reached. This being so, they manifestly 
form a compressible system. In this experiment, the information dimension Di, and the Lyapunov 
exponents of the coagulated floaters, is measured. The trajectories and the velocity fields of the 
particles are captured in a sequence of rapidly acquired images. Then the velocity sequence is 
randomly shuffled in time to generate new trajectories. This analysis mimics the Kraichnan ensemble 
and yields properties of a velocity correlation function that is delta-correlated in time (but not in 
space). The measurements are compared with theoretical expectations and with simulations of 
Boffetta et al., that closely mimic the laboratory experiment reported here. 

PACS numbers: 47.27.-i, 47.27.ed, 47.52. +j 



INTRODUCTION 

It has become a standard procedure to study fluid 
flows by tracking the trajectories of small particles intro- 
duced into the flow. Assume first that the fluid is highly 
viscous and moves azimuthally between two concentric 
cylinders, one of which is rotating at constant angular 
velocity. A series of consecutive and rapid velocity mea- 
surements will reveal that each particle is moving in a 
circle. Now imagine that the series of consecutive veloc- 
ity measurements is randomly rearranged, as one might 
shuffle a deck of cards. If the shuffled velocity field is then 
used to evolve tracer particles (explained in the Exper- 
iment section), they will appear to be jerking back and 
forth. Nonetheless, the circularity of the trajectories will 
be recognizable. Now imagine that the fluid is turbulent 
rather than laminar and that the frame rate of the cam- 
era (in Hz) is much faster than the inverse lifetime of the 
smallest eddies. Will the randomly shuffled sequence of 
measurements contain information about the turbulence, 
as in the laminar case? 

This laboratory study was motivated by a computer 
simulation of this same problem by Boffetta et al. (to 
be called BDES) Q. That simulation, in turn, was in- 
spired by a theoretical analysis of velocity fluctuations by 
R. Kraichnan [2| . Concentrating on the two-point veloc- 
ity autocorrelation function, Kraichnan argued that one 
could increase our understanding of turbulence by mak- 
ing the simplifying assumption that velocity fluctuations 
are delta-correlated in time, while retaining the spatial 
dependence of the fluctuations. Kraichnan's approach is 
related to prior work of L. F. Richardson, who analyzed 
the time evolution of the separation of particle pairs in 
turbulent flows by treating their motion as Brownian-like 

This experiment, and the computer simulations of 
BDES, focuses on an unusual type of flow, namely the 
motion of low-density particles that float on a turbulent, 
three-dimensional incompressible fluid (water in these ex- 



periments). Even though the underlying fluid is incom- 
pressible, the floating particles form a compressible sys- 
tem. If the water molecules should move downward at 
some point at the surface (z = 0), the floaters, which can- 
not follow this downward motion, will accumulate there. 
Likewise they will flee the upwelling points. Since the un- 
derlying flow is assumed to be turbulent, the pattern the 
floaters form will fluctuate in both time and space. The 
floaters are chosen to be small enough that they follow 
the velocity field at the surface; their inertia is negligi- 
ble (see below). Neglecting the wave motion of the sur- 
face (which is negligible), the floaters move in a plane. 
They are merely sampling the velocity field of the in- 
compressible flow at the surface, and hence their velocity 
at each instant t is the same as that of the surface wa- 
ter molecules (v x (x, y,0, t),v y (x, y, 0,t)) = v(x, i)). The 
geometrical pattern formed by the floaters has a fractal 
appearance, as seen in Fig. [2] (a), which reaches a statis- 
tical steady state in approximately Is. This work focuses 
on the geometric and dynamic behavior of the floaters in 
the steady state. 

The coagulation phenomenon studied here and in 
BDES has nothing to do with the interaction between 
the floaters or the waves that are present on the surface; 
all the effects observed in the laboratory are also seen in 
computer simulations, where these two effects are not in- 
cluded [H, HI]. If the flow field is isotropic, one can define 
a dimensionless compressibility at the surface, 
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which lies between (incompressible flow) and 1( poten- 
tial flow). In the simulations C is an adjustable param- 
eter, but in the experiments it is measured to be very 
close to 1/2 jij. The computer simulations yield a value 
very close to this [H, HJ . In these experiments a descrip- 
tor of the coagulation pattern is extracted from the data, 
namely the "information dimension" D\ , which is defined 
below. It is simply related to C in the Kraichnan model 
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Thus, the measurement of D\ provides an experi- 
mental test of the applicability of the Kraichnan model 
for temporally decorrelated data. 

EXPERIMENT 

The experiments are carried out in a tank lm x lm 
in lateral dimensions, filled with water to a depth of 
30cm. Turbulence is generated by a pump (8hp) which 
re-circulates water through a system of 36 rotating jets 
placed horizontally across the tank floor (see Fig. [TJ). 
This system for generating the turbulence ensures that 
the source of turbulent injection is far removed from the 
free surface where the measurements are made More 
importantly, the method also minimizes the amplitude of 
surface waves, which are unavoidable. Prior experiments 
have demonstrated that the waves at the surface have 
negligible influence on the experimental results |4j . 

The hydrophilic particles chosen here are subject to 
capillary forces which are very small compared to forces 
coming from the turbulence, and do not affect the results 
as they do in 5, 6]. The non-inertial character of the 
particles is minimal because the product r s Ai, where Ai 
is the largest Lyapunov exponent and t s is the stopping 
time of the particle 0j, is of the order 5 x 10~ 4 , which is 
too small to lead to any significant inertial effect 

During an experimental run, floating particles (50/xto 
diameter and specific gravity of 0.25) are constantly 
seeded into the fluid from the tank floor, where they 
undergo turbulent mixing as they rise due to buoyancy. 
Their motion is tracked via a high-speed camera (Phan- 
tom v.5) situated above the tank. The camera field-of- 
view is a square area of side length L — 9cm. The con- 
stant particle injection is necessary to replace particles 
at the surface during the experiment. The source and 
sink structure at the surface fluctuates in both time and 
space, which can cause particles to leave the camera's 
field of view. This experimental scheme has been used 
several times before [J H| , including a detailed compari- 
son of the experiment with numerical simulations. 

The instantaneous velocity field of the floaters is mea- 
sured using an in-house developed particle imaging ve- 
locimetry (PIV) program, which processes the recorded 
images of the surface particles. The constant injection of 
particles ensures that surface sources and sinks receive 
an adequate coverage on the surface. The local particle 
density at the surface determines the average spacing of 
the velocity vector fields produced by the PIV program. 
The resulting velocity vectors are spaced (on average) by 
Sx = 2.5 r\ over both sources and sinks, where r\ is the 
dissipative scale at the surface (Table I). This density of 
the vector grid spacing is important for the Lagrangian 
particle evolution scheme, which is discussed below. 

The camera height above the water surface was chosen 
so that a pixel size is roughly 0.1 mm, comparable to 



r\ (see Table I). Data were taken for several values of 
Re\ ~ 150 — 170 (defined in Table I) with an average 
Re\ ~ 160. The measurements in this experiment show 
no systematic variation with the Rex over this range, so 
they are used as an ensemble average so as to decrease 
measurement errors. Turbulent parameters measured at 
the surface are listed in Table HI The camera frame rate 
is set between 100-200 Hz, such that each instantaneous 
velocity vector field is separated by approximately 4— 8r r) , 
where t v is the lifetime of eddies of size r\ (see Table H . 
Results reported here do not depend on this frame rate. 

The experimentally measured velocity fields were then 
used to solve the equation of motion for Lagrangian par- 
ticles : 

^=v(x^),i), (2) 

where v(xi,i) is the velocity field and Xi = (xi,yi) are 
the individual particle positions. 

To achieve accurate results for the Lagrangian particle 
evolution, the vector fields used in Eq. [5] were interpo- 
lated from the experimentally determined velocity vec- 
tors via a bi-cubic interpolation scheme developed for nu- 
merical simulations, as discussed in Q and implemented 
in [4]. This scheme uses the smooth flow between grid 
points separated by length scales comparable to r\ to in- 
terpolate the velocity field between measured grid points. 
To use this scheme it is necessary for the measured grid 
spacing to satisfy the criterion Sx < ■nr\, where Sx is the 
measured velocity grid spacing. We have tested to en- 
sure that the results do not depend on the velocity grid 
spacing by varying the spacing from Sx — 2.5 77 to Sx = 4 
r\ (increasing this average spacing even further degrades 
the statistical quality of the data) . 

For the correlated flow data, a uniform distribution of 
Lagrangian tracers is generated at time t/ro = 0, where 
t = SU(lo)/lo is the large eddy turnover time. Here, 
To = 0.4s. The uniform distribution is then evolved using 
the above-described scheme with the velocity vector field 
record in the order it was taken during the experiment. 
To de-correlate the data, we take successive snapshots 
of the Eulerian velocity field and put them in random 
order. An initially uniform particle distribution is then 
evolved from this shuffled record. The velocity vector 
field record spans roughly 60 tq . The results in this work 
were unaffected when this record was cut in half, suggest- 
ing that the data record is long enough to mimic delta 
correlation in time. Visualizations of the correlated and 
decorrelated particle distributions (in the steady state) 
are shown in Fig. [2fa) and (b) respectively. The qualita- 
tive impression from the two images is that both distribu- 
tions lie along line-like structures, but for the correlated 
flow these structures span up to several integral length 
scales, lo/rj ~ 70, where Iq is defined in Table I. 
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FIG. 1: Schematic of the top-view (top panel) and side-view 
(bottom panel) of the experimental setup. 36 rotating capped 
jets are placed horizontally on the tank floor (shown as ran- 
domly oriented Z-shaped patterns) that pump water into the 
tank re-circulated by a 8hp pump. The central region of the 
water surface is illuminated by a laser-sheet. A high-speed 
digital camera suspended vertically above this central region 
captures images of the light scattered by buoyant particles 
(50 fim hollow-glass spheres of specific gravity 0.25). 



INFORMATION DIMENSION Z>i 

The information dimension D\ is measured by dividing 
the system's field of view into boxes of size r — s/rj, where 
s is the box size in cm, making r dimensionless. One then 
calculates the probability Pi(r) of box i being populated. 
The information dimension D\ is defined as: 



Di = lim 



-I(r) 



log r 



(3) 



The function /(r), is called the information function, 
is defined as 



(4) 
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FIG. 2: Visualization of particle distributions for the cor- 
related flow (a) and the decorrelated flow (b) at t/ro — 5. 
The time scale is in units of the large eddy turnover time 
to. The initial distribution (t/ro = 0) is homogeneous with 
4xl0 5 particles. Both particle distributions are in the steady 
state. Notice that for the correlated flow, there are coher- 
ent structures at large scales (r > 50). For the decorrelated 
flow, these large scale structures are not present, but line- like 
coagulations persist at scale r < 10. 



To calculate £>i, one plots I(r) vs logr. The slope of 
this line is D\ over the interval of r where this plot is a 
straight line. Fig. [3]shows I(r) vs. logr of the correlated 
and decorrelated particle distributions seen in Fig. [2] 
For the correlated flow, D\ = 1.16 ± 0.02 and for the 
decorrelated flow, D\ = 1.09 ± 0.02, measured over the 
interval 0.2 < r < 3. We choose to label this range of r 
as being in the dissipative range. This is consistent with 
the fact that the transition from the inertial to dissipative 
range in many experiments occurs at r > 1 (ToL lllj . 

The qualitative impression that both distributions lie 
along line-like structures is consistent with D\ being close 
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TABLE I: Turbulent parameters measured at the surface. 
Measurements are made at several values of the Rex with 
an average Re\ ~ 160. The parameters listed are averages, 
with deviations less than 10%. 



Taylor microscale 




0.47 

(cm) 


Taylor Re^ 


R v rms X 


160 


Integral Scale 


, _fj» ("II 0+*>|| (»)) 


1.42 


'° ((,,(«))«} 






(cm) 


Large Eddy Turnover 
Time 


T ° = 


0.43 (s) 


Kolmogorov Time 


r IJ = - 1/2 


0.04 (s) 


Dissipation Rate 




6.05 

(cm 2 /s 3 ) 


Kolmogorov scale 


*=(-) 


0.02 

(cm) 


RMS Velocity 


V rms = \J {v 2 ) — (v) 2 


3.3 

(cm/s) 


Compressibility 




0.49 ± 
0.02 
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FIG. 3: Plot of information function I(r) versus r — s/r) 
for the correlated and decorrelated flows at t/ro = 5. Scale- 
free behavior is observed for 0.2 < r < 3. The inset shows 
the quality of the scale-free behavior. For the correlated flow 
£>i = 1.16 ± 0.02, and for the decorrelated flow Di = 1.09 ± 
0.02. The most conspicuous difference between the correlated 
and decorrelated flows is the lack of scaling in the inertial 
range for the decorrelated flow. This reflects the loss of large- 
scale structures when the flow is decorrelated. The prediction 
for Di in the Kraichnan ensemble is Di ~ 1. 



to unity. This finding is consistent with the DNS per- 
formed in [H, where it was found that D\ = 1.15 for 
the correlated flow. At values of r > 10 (in the iner- 
tial range), the correlated and decorrelated distributions 
differ significantly. For the correlated flow in Fig. 2 (a), 
the line-like structures extend over several integral length 
scales Iq. These large-scale structures of the correlated 
flow give rise to a value of D\ that is larger than its 
small-scale value (0.2 < r < 3), as more clearly seen in 
the inset of Fig. 3. This inertial range scaling is consis- 
tent with the results of a model for adjustably compress- 
ible flow in 12j. For the decorrelated flow, the distribu- 
tion is more homogeneous over the inertial range, where 
scale-free behavior is not apparent. This demonstrates 
(quantitatively) the importance of time correlation for 
the formation of large-scale coherent structures. 

The Kraichnan model of time-decorrelated velocity 
fluctuations provides a simple relation between the in- 
formation dimension D\ and the compressibility C [l|: 



Di 



1 + 2C 



(5) 



both of which are measured in this study. The transi- 
tion to strongly compressibile flow takes place at C — 
1/2, which is approximately what is seen in this ex- 
periment. For the decorrelated flow in this experiment, 
Di = 1.09±0.02, which is measurably greater than unity. 
The decrease in the dimension D\ for the decorrelated 
flow in this experiment is consistent with the results in 
[lj], but the value obtained is larger than the Kraichnan 



prediction of D\ ~ 1. 



LYAPUNOV EXPONENTS 

The information dimension D\ was obtained in BDES 
from a measurement of the two Lyapunov exponents, 
which have been measured in this experiment as well. 
The Lyapunov exponents of the floaters give a measure 
of the mean rate at which particle pairs separate on the 
surface. In d — 2 dimensions, the largest is Ai (which 
is positive), and the smallest A2 < 0. Since the floaters 
coagulate, their sum is negative. These exponents are 
defined so that each pair is initially separated by a dis- 
tance ro that is very small (less than 77 in the turbu- 
lent case). Since the Lyapunov exponents define an av- 
erage rate of separation, a pair must be tracked for a 
long enough time that their rate of separation has be- 
come time- independent, but not so long that their final 
separation has become greater than 77 [131 ] . 

Instead of allowing particle pairs to separate contin- 
uously, Ai is calculated by the algorithm given in [Tij ]. 
After each frame, a given particle pair's separation r(t) 
is reset to ro and evolved again. This is done in the limit 
that t — > 00, where Ai should saturate to it's steady- 
state limiting value. The particle pairs are evolved using 
the same Lagrangian evolution scheme discussed above. 
The results presented here do not depend on ro, which 
was varied from 0.25?? to 1?/. An initial separation of 
ro = 0.25/7 is used. Here, we evolve an ensemble of parti- 
cle pairs (10 6 ) for a maximum of 5to, which corresponds 
to roughly 100r ?; (to =0.4 s and t v — 0.04 s). The calcu- 
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TABLE II: Parameters measured in this experiment, the sim- 
ulation of BDES, and a prior experiment where (dS/dt) was 
measured [31 • 





Di 


Ai 


A 2 


(dS/dt) 


Experiment 




Correlated 


1.16 ± 
0.02 


0.6 ± 
0.05 


-3.75 ± 
0.26 


-3.15 ± 
0.5 


Decorrelated 


1.09 ± 
0.02 


0.5 ± 
0.05 


-5.2 ± 
0.44 


-4.7 ± 
0.8 


BDES Simulation 




Correlated 


1.15 








Decorrelated 


1.05 








Bandi Experiment 








-2.4 ± 
0.02 



lation of Ai is thus: 



A 1= lim i(log^), 



(6) 



where the angular brackets denote an ensemble average 
over the various particle pairs. From Fig. |4] and [5] it is 
seen that after a small transient time (of roughly To), Ai 
saturates to a value of 0.6 Hz (correlated) and 0.5 Hz 
(decorrelated) (see Table HI)) . Also shown in Fig. @]and 
[5] are values of D\ measured at various times during the 
experiment after Ai has saturated. 

One can calculate the value of A2, using Ai and D±, 
from the Kaplan- Yorke relation (which is exact in two 
dimensions) [15j : 



£>i = 1 + Ai/|A 2 



(7) 



From this equation, A2 is calculated for both the cor- 
related and decorrelated flows in Fig. [4] and El respec- 
tively. The results from this experiment and from the 
simulations in BDES, which are summarized in Table HJ 
confirm what the eye sees in Fig. [2j D\ is roughly equal 
to unity (the particles cluster into string-like clusters), 
and the sum of the Lyapunov exponents is negative. 

There is an additional reason for measuring Ai and A2; 
their sum is related to a quantity measured in a prior 
experiment, providing an additional check on the inter- 
nal consistency of all the experimental observations. In 
a prior experiment by Bandi et al. Q, the mean rate of 
change of the entropy of the floaters, (dS/dt), was mea- 
sured in the same tank used in the present experiments, 
though at a Re\ of 100 instead of 160. Here S is the 
entropy, defined as S = — J2i Pi m (-Pi)> which is Eq. (|4]) 
in the limit that the box size r is sufficiently small. 

This mean entropy rate (dS/dt) is the spatially- 
averaged divergence of the velocity field at the surface 
[§]. It is also equal to the sum of the Lyaunov exponents, 
(dS/dt) = A1+A2, if the turbulent flow is statistically sta- 
tionary [l^]. In Q (dS/dt) ~ —2.4 Hz, while the present 
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FIG. 4: Calculation of the largest Lyapunov exponent Ai for 
the correlated flow. The information dimension Di is cal- 
culated (inset) at several instants of time once the largest 
Lyapunov exponent has saturated to a steady-state value of 
the Ai = 0.6Hz. From Eq. one can estimate the value of 
the second Lyapunov exponent A2 = —3.5 ± 0.26Hz from Ai 
and Di . The sum of the Lyapunov exponents is related to the 
entropy rate S — X1+X2 = —2.9 ±0.3Hz. This is comparable 
to the entropy rate S = —2.4 ± 0.02Hz measured in [8j. 
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FIG. 5: Calculation of the largest Lyapunov exponent Ai for 
the decorrelated flow. The information dimension D\ is cal- 
culated (inset) at several instants of time once the largest 
Lyapunov exponent has saturated to a steady-state value of 
the Ai = 0.5Hz. From Eq. one can estimate the value of 
the second Lyapunov exponent A2 = —5.2 ± 0.44Hz from Ai 
and Di . The important difference between the correlated and 
decorrelated flows is not the change in the Lyapunov expo- 
nents, but the decrease in the dimension D\ , which is consis- 
tent with the results in [J 



experiment gives E = —3.15 ± 0.5 Hz. In light of the 
uncertainty in the measurement of the entropy rate in 
@, and the difference in Re\, this agreement between 
the two types of measurement of S seems quite satisfac- 
tory. In comparing the simulation of BDES with these 
experiments, it should be born in mind that the labora- 
tory results have appreciable uncertainty, while there is 
no error (or units) quoted for the simulations. 
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SUMMARY 

Turbulence is especially difficult to understand be- 
cause the velocity fluctuations, which lie at its heart, are 
chaotic in both space and time, with very many degrees 
of freedom being excited. The phenomenon is especially 
intractable when viewed in the Eulerian frame, where 
the small-scale velocity fluctuations are contained in, and 
moved about, by eddies of larger size. 

In an effort to capture the essence of turbulence, R. 
Kraichnan, and L. F. Richardson before him, proposed 
that the problem is more easily understood in the La- 
grangian frame, where large-scale sweeping effects are 
absent. To effect further simplification, Kraichnan pro- 
posed that the temporal velocity fluctuations be treated 
as in Brownian motion, i.e., that the correlation time 
of eddies of all sizes, be approximated as being infinitely 
short (a Gaussian approximation, as in Brownian motion, 
is made as well). 

This experiment probes a type of compressible flow, 
namely the chaotic motion of particles that float on 
an incompressible and strongly turbulent fluid. The 
spatio-temporal dynamics of these floaters is studied, and 
then the observed velocity fluctuations are randomized in 
time, so as to mimic the simplifying dynamics analyzed 
by Kraichnan. As predicted, certain aspects of the par- 
ticle motion and topology are retained in the temporally 
decorrelated field. 

The topology of the decorrelated floaters is string-like, 
just as in the correlated case (compare the two panels 
in Fig. [2]) The information dimension D\ remains close 
to unity, at least at small spatial scales where this pa- 
rameter is defined. On large scales, temporal decorrela- 
tion truncates the length of the string-like structures, as 
seen in Fig. [5J The present measurements, and those of 
BDES 1], show that D\ (and presumably other fractal 
dimensions) is measurably altered by the temporal ran- 
domization, but at small scales the particle distribution 
retains its string-like character. 

The effect of temporal randomization on the two 
Lyaponuv exponents is more delicate. For both the cor- 
related and decorrelated flows, the positive exponent Ai 
remains positive (patches of particles are stretched into 
strings in both cases), but decorrelation compacts the 
structures into thinner strings, i.e., A2 is decreases from 
roughly -3 to -5 Hz. 



The steady-state entropy rate of change dS/dt, which 
for a flow admitting SRB statistics is the sum of Ai and 
A2 (la, LL7J , is lowered by decorrelation; it becomes even 



more negative. The measurements reported here, and the 
simulations of Boffetta et al., suggest that our under- 
standing of turbulence can be enhanced by comparing 
and contrasting correlated and temporally decorrelated 
flows. 
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